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Abstract 

We analyze the pattern of normal modes in linearized Lorentz- violating massive gravity over the 5-dimensional 
moduli space of mass terms. Ghost-free theories arise at bifurcation points when the ghosts get out of the 
spectrum of propagating particles due to vanishing of the coefficient in front of u> 2 in the propagator. Sim- 
ilarly, the van Dam-Veltman-Zakharov (DVZ) discontinuities in the Newton law arise at another type of 
bifurcations, when the coefficient vanishes in front of k 2 . When the Lorentz invariance is broken, these 
two kinds of bifurcations get independent and one can easily find a ghost-free model without the DVZ 
discontinuity in the moduli space, at least, in the quadratic (linearized) approximation. 



1. Introduction. The theory of massive gravity [TJ [5] attracts a new attention these days [31 SI because 
of the growing belief in acceleration of Universe expansion (the "dark energy" phenomenon) [5]j^] However, 
the violation of general covariance in massive gravity is long known to produce a number of non-trivial effects 
like occurrence of ghosts and the lack of perturbative regime at small distances [8], moreover, the van Dam- 
Vcltman-Zakharov (DVZ) discontinuities [9j [10] and the Boulware-Deser instabilities [IT] arise whenever one 
tries to eliminate the ghosts. In fact, these problems can probably be avoided, if one sacrifices the Lorentz 
invariance |12l [3J , what allows to extend the number of possible mass terms and go around the most unpleasant 
singularities in the moduli space. This was demonstrated at the level of the linearized gravity with quadratic 
actioiJ3 

(k^karifjv + k^kprjav + Kkarip^ + k u kpr/ ai _^j - (k^riap + kakpi]^ - ^(r/^rivp + ^7]^ + fc 2 ?^/^} Wh' 

+ mlhl + 2m 2 1 h 2 Qi - m\h^ + m^h^ - 2m\hQ h li , (1) 

where the first line is nothing but quadratic approximation to the Einstein-Hilbert action, while the second line 
contains five different mass termsj^] which violate both gauge (general coordinate) and Lorentz SO(d — 1,1) 

invariance, but preserve space rotation symmetry SO(d — 1). In our notation, = ^) > while 

h 2 j — Y^,i hfj- The theory also has the P and T reflection symmetries, so that all scalar physical quantities 
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1 As we understand, the argument is as follows. There is an experimental evidence that cosmological constant is actually non- 
vanishing. From the point of view of flat geometry, the cosmological constant makes graviton massive (in fact it also provides it with 
a source term, linear in h, therefore a more accurate analysis, including the change of expansion background, is actually required), 
then, Lorentz-invariant massive gravity looks ill, but this can be cured by switching on especially adjusted Lorenz-violating terms, 
perhaps as small as the cosmological constant, what makes their effect small and consistent with existing observations. Alternatively 
one could say that non-vanishing negative cosmological constant implies that flat background geometry is substituted with AdS 
one, and all analysis should be made differently from this new perspective [6]. In fact motivations for the study of infrared-modified 
gravity are not exhausted by the dark-energy problem, for some other examples see [3] and [7], 

2 Throughout the paper, our convention for the metric signature is (—,+,..., +). 

3 Of course, one can violate Lorentz invariance not only in the sector of masses, but also in kinetic term and add higher 
derivatives in space directions, which do not produce new ghosts. For profound example of this kind see [13] ■ The methods of the 
present paper are straightforwardly applicable to these non-minimal deformations, but on this road the moduli space M is in no 
way restricted and eigenvalue patterns can be made arbitrarily complicated. 



depend on the squares uj 2 and k of frequencies and space momenta. The Lorentz invariance is restored if 
the five mass parameters can be expressed through only two independent quantities, A and B: 

ml = B — A, 

,2 _ „,2 



m\ = m\ = B (2) 



m 1 = m 2 = A, 

and ZC M „ iQ /3 in (|T|) reduces to 

— {kfikaTjp^ + k^kprj av + k v k a r]py, + k v kpf\ al ^ — {^k^k u rj a f} + kakpr)^ 

-^(fc 2 + A) (jl^aTlvp + llvaVnf^) + {k 2 + B)rj^ v T] a p (3) 

The ghost-free Lorentz-invariant Pauli-Fierz [T] massive gravity corresponds to the choice A — B: it, 
however, suffers from all the above-mentioned problems and thus looks unviable [3]. The Lorentz- violating 
theories ([lj can be ghost free when either mo = or mi = 0, and the second choice is the current favorite 
candidate for a phenomcnologically acceptable version of massive gravity [3J. 

Lorentz violation breaks a lot of familiar properties of quantum field theory models and looks unusual in 
many respects. It gives rise to the whole variety of non-trivial quasi-particles which can be ghosts, superluminals 
and even not look like particles at all (cither relativistic, or non-relativistic) . In |14j we provided a systematic 
analysis of the theory and carefully reproduced and explained the results of [3J , also relating them to the 
obvious self-consistency of Kaluza-Klein theories, which involve massive gravitons but remain free of any kind 
of pathologies. Here we present this analysis in still another, concise and formal way, omitting a lot of details 
and physical motivations included into [U]. Note that, due to the different choice of signature, the sings of 
eigenvalues throughout the paper are opposite to those in [T4"] . 

2. The main quantity: the propagator II(fc) over the moduli space. We remind briefly the stan- 
dard string-theory approach to consideration of a family of physical theories |15j . adapting it to a particular 
application to linearized gravity, perhaps, Lorentz- violating. 

The physical content of a particular theory (model) is best expressed in terms of the partition function 

Z(J) = J Deb e t(SW+ f ^ (4) 

In quadratic approximation, when 

S(4>) = J d d k <j){-k)K{k)<j){k) (5) 

and 

' d d x J(j>= I d d k J(k)cj)(k), (6) 



this Z{ J) is also a quadratic exponential, 

Z(J) = exp f-i J d d k J{-k)K- 1 {k)J(k)^j (7) 

made from inverse of kinetic matrix K{k), i.e. the propagator. This is a finite-dimensional matrix in the space 
of fields 4>(x): if (j> a (x) carries an index a, then K a ^(k) carries two indices a, b. In the case of vector fields a is 
just the Lorentz index fi, in the case of gravitational field a = (/iz^) is a symmetric pair of the Lorentz indices, 
thus taking d ^ d ^ different values, what reduces to ( d ~ 1 H d + 2 ) m the case of the traceless field. Our main task 
is to investigate the quantity II(fc) = J(—k)K^ 1 (k)J(k). 

Of most interest for us in this paper are two kind of characteristics of II. 

(i) A singularity of II(fc) defines a propagating particle and the position of the singularity defines its dispersion 
relation u> = s(\k\). 

(ii) The quantity V(|fc|) = II (a; = 0, k) defines an instantaneous Newton/ Coulomb/ Yukawa- like interaction. 
The partition function Z and its quadratic approximation exp(— | J II) are of course defined over the space 

of theories M. (and are, hence, generalized r-functions |16j . ordinary and quasiclassical respectively), and we 
are going to study the singularities (reshufflings or bifurcations) of dispersion relations and potentials V over 
the moduli space M.. In the current problem coordinates in M. parameterize the kinetic matrix K(k), actually, 
the mass terms, and, as usual in string theory, in the spirit of third-quantization, they can be considered as 
vacuum averages of some other fields (slow variables or moduli). 
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3. The notion of eigenvalues and its ambiguity. The problem of dispersion relations is basically that 
of the eigenvalues of K(k): roughly, id = e(\k\) is a condition that some eigenvalue A(fc) = 0. However, this 
"obvious" statement requires a more accurate formulation. The point is that K is actually a quadratic form, 
not an operator, what means that it can always be brought to the canonical form with only ±1 and at 
diagonal, thus leaving no room to quantities like A(fc). Still, this "equally obvious" counter-statement is also 
partly misleading, because we are interested not in an isolated quadratic form, but in a family of those, defined 
over M.. This means that the sets of ±1 and can change as we move along M., and degeneracy degree of 
quadratic form K(k) can change. Of course, this degree (a number of 0's at diagonal) is an integer and changes 
abruptly - and thus is not a very nice quantity. A desire to make it smooth brings us back a concept of A(fc). 
However, in order to introduce A(fc) one needs an additional structure, for example, a metric in the space of 
fields. 

In application to our needs one can introduce "eigenvalues" X(k) as follows: consider instead of n = J -4 J 
a more general quantity 

n(A|fc) = Jj J- Jj j ( 8 ) 

Then as a function of A it can be represented as a sum of contributions of different poles: 

n(A|fc) = £ (9) 

then X a (k) are exactly the "eigenvalues" that we are interested in, and our original 

v ^ a b a c (k)J b (~k)J c (k) 

n(fc) = \ (10) 

a,b,c 

The only thing that one should keep in mind is that this decomposition depends on the choice of additional 
matrix (metric) /, which can be chosen in different ways, in particular, its normalization can in principle depend 
on the point of M. . We shall actually assume that it does not, and clearly the physical properties do not depend 
on this choice, however, concrete expressions for A a (fc) do. It is important, that the dispersion relations - the 
zeroes of X a (k) - are independent of I. 

Introduction of / is also important from another point of view. To be well-defined, the Lorentzian partition 
function requires a distinction between the retarded and advanced correlators (Green functions), which is usually 
introduced by adding an infinitesimal imaginary term to the kinetic matrix K: the celebrated ie in the Feynman 
propagator]^] However, in the case of kinetic matrix this is not just ie, it is rather ielp with some particular 
matrix Ip. If we identify our / with Ip, then the dispersion relations are actually 

\ a (k) = ie (11) 

what implies that A a (fc) is, in fact, very different from — X a (k), and this is related to the important concept of 
ghosts. 

The most natural choices of the matrix / are probably either just the unit matrix, or "the Lorentzian 
unit matrix", i.e. that with -1 corresponding to the 0-components. The physically justified choice is the unit 
(Euclidean) matrix, while technically it is often simpler to work with the Lorentzian unit matrix, especially 
when dealing with theories with the Lorentz invariance unbroken. At the same time, in these two cases it is 
only ghost content of the non-scalar sectors which differs. Therefore, it is often safe (and technically preferable) 
to use the Lorentzian unit matrix. We illustrate this in the simplest warm-up example of the massive vector 
field theory in Appendix I, where we compare the results obtained for the two cases of Euclidean and Lorentzian 
eigenvalues (Euclidean and Lorentzian unit matrices). Since this paper is rather devoted to the method than 
to concrete physical applications, we use the Lorentzian eigenvalues here, leaving the Euclidean ones for 14J, 
where we deal with physical issues. 



4 The Feynman propagator implies that particles with the dispersion relation u = +e(|fc|) propagate forward in time, while 
antiparticles with u> = — e([fc|) - backwards in time. Since 9(±i) = J ± u -^ie i we have for the propagator 

1 / 1 1 \ _ 1 

2e \u! — e — ie —u) — e — ie J u) 2 — e 2 — ie 



For ghosts with the propagator ^ 2 _* 2+ie the situation is inverse: particles propagate backwards while antiparticles forward in 
time. See also Appendix II. 
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4. Spectrum and the phase diagram. Important information about the theory is contained in its spectrum: 
positions of the poles of II(A|fc) in the complex A-planc. These positions define the dispersion relations A a (a->, k) = 
between the frequency w and the wave vector k of elementary excitations (quasiparticles) and the way these 
relations depend on the point of the moduli space M. . 

As one knows well from condensed matter physics, in generic Lorentz-violating theory dispersion relations 
are quite sophisticated, they are roots of polynomial equation and often do not possess any useful analytical 
expressions. Sometime they are better represented by pictures: the plots X a (Lo) or A a (fc), however when the 
pattern is multi-dimensional one can only draw its particular 2d or 3d sections, which do not provide complete 
visualization. Developed algebra-geometric intuition is actually needed to analyze the spectrum - surprisingly 
enough, this is already the case in such a fundamental (and seemingly simple) theory as linearized gravity! 

Of main interest are qualitative features of the spectrum and their bifurcations: the changes of these qual- 
itative features when one goes from one region of the moduli space to another. The corresponding division 
of the moduli space into domains with qualitatively different spectra (and, perhaps, other physically relevant 
characteristics like structure functions a^ c (k)) is called the phase diagram of the theory (or, better, of the family 
of theories) . 



5. Ghosts, tachyons, superluminals and DVZ jumps The simplest examples of qualitative features 
of the spectrum are the presence or absence of exotic (from the perspective of Lorentz-invariant field theory) 
excitations, like ghosts or superluminals. 

The ghost differs from the normal particle by a sign in front of lu 2 in X a (k). For example, for a scalar particle, 



d\ a (k) 



du> 2 



A a (fe)=0 



< normal particle 
> ghost 



(12) 



In order to define this sign one needs to compare it with the one in front of ie in (111, see Appendix II for 
a more detailed discussion. Problems are actually expected when excitations with opposite signs are present: 
when the "ghosts" need to coexist and interact with the "normal" particles. The condition 



d\ a {k) 



du 2 







(13) 



A a (/c)=0 



defines the loci in the moduli space, where the ghost content of theory can change. 
A similarly-looking condition 



9A a (w = 0,fc) 



dk 2 



(14) 



A o (fc)=0 



defines the loci of the DVZ jumps, see below. In the Lorentz invariant theory, where A a depends on k 2 — 
—ijj 2 + k 2 , the two conditions (12) and (14 1 are clearly related. Therefore, one can easily come across the DVZ 
jump when trying to get rid of ghosts - and this, indeed, happens in the simplest Pauli-Fierz version of linearized 
gravity. After the Lorentz violation, the link between (12) and (14 1 is relaxed. 



Next, the difference between the normal particles and tachyons is as follows: 

if A a (w, k = 0) = has real solutions for the frequency w, this is a normal particle, 
if \ a {oj — 0, k) — has real solutions for the wave vector fc, this is a tachyon. 

The superluminal propagation [17] is controlled by the group velocity 

_ d\ a {k)/dk 



d\ a (k)/du 



(15) 



(16) 



A a (fc)=0 



in the usual way: 



< l 
= l 
> l 



normal particle 
light — like particle 
superluminal particle 



(17) 



In fact, it makes sense to further distinguish between different superluminals by looking at another quantity 

d\ a (k)/dk 2 



v 2 = 



d\ a {k)/duo 2 



(18) 



Aa(fe)=0 
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For the ordinary relativistic particle with A = — lo 2 + k 2 



this V 2 = 1 independently of the value and even of 



the sign of mass m 2 . Thus, the ordinary tachyons with negative m 2 and v 2 > 1 are rather "soft" superluminals. 
In the Lorentz violating theories things are much worse: there are "harder" superluminals with V 2 > 1. 

Finally, the DVZ jump can occur when one of the scalars becomes infinitely heavy. Then the massless 
limit, when all the five moduli Too, . . . , to 4 — > 0, gets ambiguous: this scalar can either remain infinitely heavy 
or acquire a finite mass or become massless, depending on a particular way the limit is taken. Thus, the 
contribution of such a scalar to the instantaneous potential is also ambiguous and depends on the way one 
approaches the point too, . . . ,TO4 = 0: if we are interested in physically relevant quantities, this point in M. 
is, in fact, singular and should be blown up to resolve the singularity. For the generic dispersion relation the 
role of mass in above reasoning is played by the root k 2 of the equation X a (oj = 0, k 2 ) = (the real mass gap 
arises when the root is negative, k$ < 0). The DVZ ju mp can occu r w hen k 2 . — > — oo, and this actually requires 
that X a (uJ = 0,k 2 ) has an asymptote which satisfies (14l. Thus, (14| is a necessary condition for a the DVZ 
jump to occur. Note, however, that (14 1 is more restrictive, because uj = condition is additionally imposed: 



thus it defines a codimension-one subspace in the moduli space M., while (12 1 can hold for particular u and 
k in codimension-zero domains of M.. The DVZ jump is basically a non-commutativity of the limits, i.e. the 
difference between the two naive definitions of the static potential (the instantaneous Newton/ Coulomb/ Yukawa 
interactions) at a given point Mq in the moduli space. Such a difference can occur when the number of degrees 
of freedom changes at Mo, i.e. when the two bra nch es of dispersion relations merge or intersect. This happens 
if the two roots A a (w = 0, k) coincide, i.e. when ( 14 1 takes place. 



In Appendix I we thoroughly study conditions of emergency of ghosts, tachyons, superluminals and the DVZ 
jumps in the example of massive vector theory. 



6. Eigenvalues and discriminant analysis. The "eigenvalues" X a (k) are roots of the characteristic poly- 
nomial 



dog C 



C/(A) = discriminant^ (s(<f>) - \((/>I(f>j) = det(K - IX) = J[ (A - A a ), 



(19) 



since the discriminant of a quadratic form is actually a determinant of the corresponding matrix. On-shcll 
conditions A Q (fc) = are zeroes of discriminant ^ (^S((j))j = dct K itself and do not depend on the choice of I, as 
we already mentioned. 



Similarly, conditions like (13 1 and (14 1 are zeroes of the ratio 



resultant 



x (c(X),SC(X) 



resultant 



(C(X),C'(X) 



(-) 



1+dcg c 



Jj6X a 



(20) 



where S is any variation of the coefficients of C(A), say, resulting from an infinitesimal change of uj 2 or k 2 , 
and C"(A) is A-derivative of C(A). Note that the resultant in the numerator has degree deg(C) = #(a) in 
the coefficients of SC, and the resultant in denominator is actually a discriminant of C(A). For definitions of 
resultants and discriminants see, e.g., [18| [19j. 



7. The pattern of eigenmodes for Lorentz-violating gravity. After these general remarks, we return 
to the concrete model: the linearized massive gravity (|T|). 

Eigenvectors of the kinetic matrix are naturally split into three groups: traceless tensors, vectors and scalars. 



In more detail, the 



d(d+l) 



components of symmetric tensor hn U are decomposed as follows: 



d(d+l) 



(d-2)(d+l) 



massive spin 2 



(d-1) + 1 + 1 

space-time transverse secondary space-time trace 
Stucckclbcrg vector 



(21) 



d(d - 3) 



+ (d - 2) 



v spatial— transverse tensor 



spatial trace 



d-2 



spatial — tr 
Dtucckclbci 



longitudinal 
icckclbcrg scala 



Stucckclbcrg 



where the first line is the SO(d — 1) classification in the rest frame (where k — 0), while the second line is the 
classification in the arbitrary frame, i.e. that w.r.t. SO(d— 2), which acts in the hyperplane transverse to k. 
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Accordingly the characteristic polynomial C(A) in the generic frame is decomposed as 

4 

C(A) = (A - A gr )^P 2 (A) d - 2 Q 4 (A) = (A - A gr )^(A - A+ c ) d " 2 (A - A"^" 2 J[ (A - A° c ) (22) 

a=l 

where P2 and Q4 are polynomials of degree 2 and 4 respectively and all their coefficients as well as X gr are 
quadratic functions of u> and fc. Since the coefficients of C are quadratic functions of lu and fc, this means that 
X gr is some bilinear combination of uj and fc, 

Aj ec = P2 ± v^4j where P2 and ^4 are respectively quadratic and quartic in lu and fc, 
Ag^ 2 ' 3 ' 4 are the roots of degree-four polynomial. 



In the rest frame the roots should be grouped in a different way, according to the first line in (21 1 


C RF (X) = (A - X gr ) ( 2 + ' (A - X vec ) d \X - A+ )(A - A sc ) 


= (23) 

fc=0 


i.e. at fc = 




A+ c (fc = 0) = A 9r (fc = 0), 
X s s f{k = 0) = A gr (fc = 0), 
X s s f (fc = 0) = A" c (fe = 0) 


(24) 


where "spT" and "spS" label the spatial trace ha and the spatial Stueckelberg scalar h Qi — fc^s respectively. The 
remaining two scalars, the space-time trace (stT) h 1 ^ and the secondary Stueckelberg scalar (seS) = k^k v a 
have eigenvalues, which are roots of quadratic equation: 


A± =q 2 ±,M\k= 


(25) 


In Lorentz invariant theory one can obtain eigenvectors and eigenvalues in an arbitrary frame by a Lorentz 
boost, but Lorentz violation forbids such a simple procedure. 


In gauge invariant theory all the d Stueckelberg fields have vanishing eigenvalues and one gets 


C G i{X) = X d (X X gr ) ( 2 ' (A - A+ c ) d 2 (A - X s s f)(X - Xf c T ) (26) 


i.e. in this case the two trace (spT and spS) eigenvalues are the roots of quadratic equation, 


A* T — t 2 ± Vu\ GI 


(27) 


Actually gauge invariant will be only the massless gravity (where, by the way, transition to the rest frame is 
not a justified operation). This can be summarized in the following scheme: 


rest frame normal modes g 


auge invariant (massless) case 


graviton — > 

/ 

massive graviton «— © 


graviton 


\ 

vector — > 


vector 


y Stueckelberg vector 
Stueckelberg (d — l)-vector \ 

\ Stueckelberg scalar — > 


Stueckelberg d-vector 


/ 

secondary Stueckelberg scalar <— secondary Stueckelberg scalar 




© «— spatial trace — > 


spatial trace 


space-time trace «— space-time trace — > 


space-time trace 
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Figure 1: The left picture plots the four eigenvalues as functions of — u) 2 in the Lorentz-non-invariant case in the rest frame and 
with 1TI4 = 0. In this case, the figure is maximally degenerated, and all the eigenvalues are straight lines. The pattern is described 
by positions of the two horizontal lines given by values of and m 2 , and by positions of the four intersections which depend also 
on m| and m|. In the middle figure, the degeneration is partly lifted by choosing non-zero 77-14 (still in the rest frame). At last, 
in the right figure, a typical perturbation of the previous figures shown, when both non-zero m? and momentum are switched on 
resolving all the four marked crossings of the left figure. The parameters here are: = —9, m 2 = — 1, m| = 9, m| = 4. 



The most interesting sector is that of scalars, with complicated inter-mixture of four eigenvectors. The 
pattern of eigenvalues is most simple at 7774 = and k = 0: 

X s t — —u 2 + m|, 
\ s = mj, 
X sS = -ml, 

\ stT = (d- 2)lu 2 + to| - (d- 1)7713 (28) 

see Figflji. 

Switching on 7714 and k leads to a bifurcation: repulsion of levels, so that Fig[lji is immediately transformed 
into Figjlj). What is important, however, the horizontal asymptotes stay at their positions: at A = — TOq and 
A = m\. 

Bifurcations become visible in the physical spectrum when one of these asymptotes coincide with the real 
axis, A = 0. Clearly, this happens when either too = or toi = 0. 

In this paper we introduce eigenvalues in a Lorentz-invariant way, taking I = II = 77 = diag(— 1, 1, . . . , 1), 
even despite the Lorentz symmetry can be explicitly violated by mass terms in the Lagrangian. This should be 
kept in mind in comparison to [14, where "Euclidean" eigenvalues were considered, associated with the choice 
I = I E = diag(l, !,...,!) (see also Appendix II). 



8. Restriction to subspace of Lorentz invariant theories in M.. To reveal the physical meaning of 
pictures like Fig(l] it is instructive to begin with the simpler Lorentz-invariant case ([2]) with only two moduli A 
and B. In this model the Lorentz symmetry expresses eigenvalues in arbitrary frame through those in the rest 
frame, so that the classification of eigenvectors is always described by the left column of the table. It remains 
only to evaluate concrete functions of k 2 = —0J 2 + k 2 : 

tensors (d+ l)(d - 2)/2 X gr = k 2 + A 



Stueckelberg vector d — 1 X vec = A (29) 

scalars 2 Xf c = A - ^ 2 +^± V{d -2n k ^ +i{d -^ 

Thus in the Lorentz invariant case we have two scalars: the space-time trace ("+" sign in above formula) and the 
secondary Stueckelberg ("— " sign). When the Lorentz invariance is violated, we will get 4 scalars, two additional 
coming from the tensor and vector multiplets (the spatial trace and spatial Stueckelberg scalar respectively). 
In Fig|2]we plot eigenvalues corresponding to these four scalars as functions of k 2 in the Lorentz invariant case 
(we include those two scalars that are parts of the tensor and vector multiplets in this case). The parameter A 
enters as a common shift of the horizontal axis, and the mass-shell condition A = for propagating particles is 
satisfied at the intersections of the eigenvalue curves in the picture with the abscissa axis. We see that the four 
lines become straight, like in Figjlji, only at B — 0, Fig|2];: this is the only point where the condition 7774 = 
is consistent with ([2]). For finite B, the two eigenvalues X sc are repulsed, like in Fig(l|D, while nothing happens 
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Figure 2: The left fi gure plots the four eigenvalues in the Lorentz-invariant case at zero momentum, and the middle figure at non- 
zero momentum (which corresponds just to shifting the whole figure to the left). The right figure corresponds to the degenerated 
case of B = 0. Note that in this case there is no way to resolve all the intersections of the right picture. 



to the other two eigenvalues, protected by the Lorentz symmetry: only two scalars can mix in this symmetric 
situation. 



We see that of these two scalars only one can be on-shell and, whatever it is, it is a ghost, see ( 12 ): the slope 
of the curve A(— ui 2 ) is negative everywhere and thus also on mass-shell, where A(— uj 2 ) = 0. It is a tachyon 
or not, depending on where the intersection with the abscissa axis occurs: to the right (tachyon) or to the 
left (normal) of the ordinate axis (for positive or negative k 2 — —uj 2 + k 2 ). The only chance for this ghost 
to disappear from the spectrum of propagating particles is when the thin line (asymptotics of the eigenvalues) 
coincides with the abscissa axis, i.e. when A = B: this is exactly the Pauli-Fierz model [1]. Clearly, at this 



point ( 13 1 is fulfilled. However, exactly at the same point in the moduli space condition (|14|) is also satisfied 



and the DVZ jump occurs (it comes with no surprise, because ( 13 ) and ( 14 1 always coincide if the Lorentz 
invariance is not violated). 

The DVZ jump occurs because the instantaneous-interaction potential V(k) — H(ui — 0, k) does not have a 
well defined limit when both A and B tend to zero. Instead, V{k) is well defined on a properly compactified 
moduli space with a blown-up singularity at A = B = 0: if one parameterize B as B = A + A 2 £, then V is 
actually a smooth function of A and £, see Fig|3}:. In more detail, the Newton/ Yukawa-potential is given by 

J (k)J (k) ( d-2 1 \ 2 2 A(dB - A) 

V{k) = (d - W - 2) l^T^ + WTW 2 ) ' m=A > = (d-2)(A-B) (30) 

see [TU s.3.6]. It is plotted as a function of A and B in Fig|3^ at some fixed value of k 2 . Poles at the two lines 
k 2 + m 2 = and k 2 + M 2 — correspond to propagating degrees of freedom, the second singularity may exist 
even at uj 2 = i.e. at positive k 2 , because the ghost can be also a tachyon, with M 2 < (it is not the case if 
B < A < dB). Clearly, the function V is discontinuous at A = B = 0, but the singularity is resolved in the 
coordinates (A : £) in Fig{3]:, at expense of gluing in a whole line (A = 0, £) instead of a single point A — B = 
(the singularity point is " blown up" ) . The situation is of course similar to the resolution of the singularity at 
A = B = in the rational function a+b passing, say, to polar or non-homogeneous coordinates B = A£', 



the difference here is that the singularity in (30 1 is rather cusp-like, see Figj3p, and the blow-up procedure is 
slightly more involved. 

9. Back to generic Lorentz violating theory. Coming back to the Lorentz-violating masses ([T]), we obtain 
a somewhat richer pattern of bifurcations, but their physical interpretations remain very similar. The essentially 
new thing is that the singular subspace of the moduli space has higher codimension and can be passed by in an 
easier way. 

We begin with the analogue of Fig|2^ in the rest frame: see Fig flp. To keep pictures similar to the Lorentz- 



invariant case, we plot dependencies on —to 2 , not +lo 2 . Instead of (|29|), we have now 
tensors 



(d+l)(d-2) 

2 "9 r 



Stueckelberg vector 



m — m 



scalars 2 Ai, = 



o-(rf-l)m3 + (d-2)w 2 ±y // (ml+ml-(d-l)ml + (d-2)u' 2 ^j + 4(d-l)r 



(31) 
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Figure 3: Resolution of singularity in potential V at A = B = 0. Plotted are potentials V at k 2 = 1 as functions of A and C, 
where C substitutes B and is introduced in three different ways: B = A + C (Fig.a), B = A + AC (Fig.b), B = A + A 2 C (Fig.c). 
Clearly, V is a smooth function at A = only in the third case. Of course, V is also singular when propagating particles contribute, 
i.e. at k' 2 + m 2 = and k 2 + M 2 = 0. These singularities is avoided since we present only a fragment of plots with small enough 
J 4<fc 2 = l. 
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-10 
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Figure 4: The eigenvalue curves at u) 2 = (which are relevant to describing the potential). The first two figures (a and b) 
correspond to the generic Lorentz violating case (the values of parameters are = —9, m 2 = —1, m| = 9, m 2 = 4, m 2 = —2 and 
rrtQ = 8, tiij = m 2 = 6, m| = m 2 = —2. The third figure (c) describes the Lorentz invariant case, when the eigenvalue asymptotics 
coincide with the abscissa axis. This corresponds to the DVZ jump and, at the same ti me, t o the Pauli-Fierz theory (A = B = 3). 
Figure d (m,Q = —9, m 2 = —3, m 2 = 9, m| = 4, m 2 = 2.4) demonstrates that condition Jl4| can be also realized in a different way 
(when one of the eigenvalue curve touches the abscissa axis). 
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Of the four crossings at PI, P2, P3, P4 in FigJTJi no one is protected by the Lorentz invariance, still in the rest 
frame only one is resolved by switching on 777,4 7^ 0, Figjljj. The remaining crossings are resolved when we also 
switch on non- vanishing fc 2 , then 



tensors dlyd 2 ^ X gr — —oj 2 + k 2 + m% 



vectors 



9y ,,_ 9l ,± _ -u J 2 + k 2 +m 2 1 +m 2 2 ±^/ (-^+fc 2 ) 2 +2(mj-m|)(^ + ^) + (m;-ml)^ (32) 

Z X (a Z) A vec — 2 



scalars 4 C4(A SC ) = 

where the polynomial C4 of degree 4 in A is explicitly presented in Appendix III, eq.(73l. The result is shown 



in FigjTJ;. In these pictures, the on-shell conditions for propagating particles correspond to intersections with 
the abscissa axis. Propagating particles disappear from the spectrum when this axis coincides with one of the 
eigenvalue asymptotics (thin lines in figures). 

In order to investigate the DVZ jumps in the instantaneous-interaction potential, one should instead look at 
FigjiJ where the four scalar A's are plotted as functions of k 2 at vanishing oj 2 . Any crossing with the abscissa 
axis in this picture corresponds to a tachyon. The DVZ jumps occur when any of the eigenvalue asymptotics 
(thin lines) coincide with the abscissa axis. For more pictures, describing the emerging phases, see Appendix 
IV. 



10. Conserved currents and instantaneous interaction. In gauge-invariant theories the currents, at- 
tached to the gauge field in Q, are conserved: if not imposed "by hands" , this condition appears automatically 
from integration over the pure gauge degrees of freedom. When the gauge invariance is explicitly broken, say, 
by the mass terms in the second line of ([l]) , this requirement is no longer enforced by the theory itself, instead it 
is imposed on massive gravity on phcnomcnological grounds: according to the currently dominating paradigm 
one is allowed to "spoil" properties of the gravity sector, but not of the matter one, which is believed to be 
under a much better experimental control. 

Conservation of currents is extremely important, because it de facto eliminates some would-be propagating 
degrees of freedom from the physically relevant quantity H(k). Let us remind that in ordinary photodynamics, 
i.e. the Maxwell theory with Lagrangian F* v , we have II = ' 2 M ' g 2 what is actually equal to 

t2 72 

n = — + J 4- (33) 

oj 2 — k 2 k 2 

for conserved current, satisfying ojJ$ = kJ = \k\ J||, so that the longitudinal photon is actually eliminated from 
LT, being substituted by a non-propagating instant Coulomb interaction. This fact persists in other theories with 
conserved currents, including ([I]), even if the gauge symmetry is violated by mass terms: the Coulomb interaction 
becomes the Yukawa one or even more complicated but continue to possess an instantaneous component. 
However, this is not explicitly seen at the level of eigenvalue analysis that we performed in this paper. This 
" drawback" can probably be cured by considering the " Euclidean eigenvalues" , suggested in [14] , but this can 
also be considered as a rather artificial trick. 

Another important remark is that even if some mode drops away from II when the currents are conserved, 
this by no means implies it can not be radiated (emitted) by a conserved current: one can easily imagine 
situations (construct models) when a mode is emitted, but can not be captured by another conserved current 
later. This happens if space-time transverse modes are mixed with the pure gauge ones - what can not be 
gcncrically forbidden in gauge- violating theories. If this happens, then the fact that the mode drops away from 
II is not sufficient to claim that it is indeed non-propagating, and one should be careful and not overlook such 
possibility. 



11. Conclusion. To conclude, we used the currently popular example of linearized massive gravity [5] to 
illustrate the general behavior of normal modes (quasiparticles) over moduli spaces of sophisticated physical 
theories and proposed to analyze this behavior by the standard techniques of linear and non-linear algebra 
[19 . Already in this relatively simple example we observe a rich pattern of bifurcations and a need to resolve 
singularities in the moduli space in order to avoid the DVZ discontinuities [H] and other pathologies. This simple 
exercise can serve as an elementary introduction to the general string theory problems from the perspective 
of ordinary - and even phenomenologically acceptable - classical field theory. At the same time, this analysis 
can help to visualize and systematize the results of [3J about the ghost-free versions of massive gravity and 
further clarify role of the Lorentz violation in constructing such a theory, at least, at the level of quadratic 
approximation. Whatever will be its relevance for phenomenological application, massive gravity looks very 
convenient for fighting prejudices of previous experience, unapplicable when the gauge and Lorentz invariances 
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are broken, and it will play a role in building new bridges between elementary particle physics and generic 
quantum field/string theory 
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Appendix I. Breaking Lorentz invariance in vector theory 

Here we consider the theory of massive vector field with the Lorentz invariance manifestly broken. 



Euclidean eigenvalues 

If one chooses the unit (Euclidean) matrix for / in d8j), one has to diagonalize the following kinetic operator: 



u 2 - (k 2 + Ml) 



K„ v = uk n uj 2 -MI (34) 



where the spatial momentum, fc|| is directed along the first direction and Mo,i are the massive terms that man- 
ifestly break the Lorentz invariance. The problem of diagonalizing this matrix leads to Euclidean eigenvalues, 
the result reads: 



A - = 2 



with 



\ (A - ^A 2 - 4(M 2 c 2 - M 2 fc 2 - M 2 M 2 )) (35) 
A+ = X - (A + .^/A 2 - 4(M 2 u 2 - M 2 fc 2 - M 2 M 2 )) (36) 



A = u? + fc 2 | + Ml - M\ (37) 



and all other d — 2 eigenvalues are equal to 

Xi = uj 2 - fcjj - Ml (38) 

There are two kinds of dispersion laws. The condition A; = evidently leads to the d — 2 excitations with the 
dispersion law 

uj 2 = fcjj + Ml (39) 
At the same time, the conditions A± = have only one solution 



^ 2 = M i+^\ (4°) 
The simplest way to see this is to look at the determinant of which is equal to 

(u? - ftjj - M 2 ^ d 2 (-MlJ 2 + M^ Ml + Mlkf\ (41) 
Now one can easily analyze these eigenvalues for physical effects: 



tachyons: Dispersion law (39 1 leads to a tachyon as soon as M 2 < 0. At the same time, dispersion law (40) 
leads to a tachyon when Mq < 0. 



(42) 



super luminal: This may come only from dispersion law ( 40 ) , which always violates Lorentz invariance unless 
Mq = Ml (since then some of the vector field modes propagate with the speed of light, and some with the 
speed of light times Mi /Mq), and, in the case of Mi /Mo > 1, describes the superluminal. 

ghosts: The ghost content of the theory is controlled by the derivatives on mass shell (i.e. at points, 
where A = 0). These are 

dh = l 
dui 2 
and 

^ = Ml (43 ) 

du 2 k 2 (M 2 +Ml) + M^ 1 ' 

This derivative is zero only when Mq. However, the ghost content of the theory can not change at this point, 
since it is Mq that enters the numerator and the ghost never emerges. If, however, Mq + M 2 < 0, there is also 
a singular point where the derivative changes the sign and, therefore, the ghost emerges. 

The condition Mq — is a counterpart of the condition niQ = in the gravity case, in this case the "live" 
excitation branch comes away from the spectrum. Another special case is Mi = when there only constant (in 
space) mode is present in the spectrum. This is an analog of the mi = condition in gravity. 



13 



DVZ jump: It is described by the derivatives Jp- at zero frequencies. The derivatives are 



Ok 2 



= -1 



(44) 



and 



d\ ± 



(45) 



Therefore, one of the derivatives is zero and, hence, there can be a DVZ jump. 



DVZ jump 



We obtained that the necessary condition of the DVZ jump requiring that (45 1 to be zero is fulfilled. However 



it is identical zero for all values of parameters which makes the general argument about the DVZ condition 
meaning lessjf] Therefore, to establish if the DVZ jump is realized, one needs a closer inspection of the interaction. 
The interaction with external currents is giv en by the term JK~ X J in the action, where J is a column ( Jq, J\, J±) 



and the propagator is the inverse of K (34 1. If one additionally requires for the currents to be conserved, 
the interaction reads 



j-2 



k\Ml 



tu 2 M 2 



k 2 kjjAf? 



uu 2 M 2 



MlM\ uj 2 



■A 



k 2 
I 



M 2 



(46) 



Bringing masses to zero in this expression in any order, as well as putting them first equal (the Lorentz-invariant 
case) and then bringing to zero leads to the same result reproducing the standard QED 



T 2 

k 2 



J2 



k 2 
II 



(47) 



If one consider a static potential in (46 1, i.e. the interaction generated by a static external current, J± 
uj = 0, one obtains 



with 



1,2 



h 2 h 2 
«H «|| 



Ml 



(48) 



which also does not shows up any jumps. Therefore, there is no the DVZ jump in this case. 

This is mostly due to a specific form of the interaction. Would be there a term, e.g., Mq instead of MqM 2 



in the denominator of (46 1, there is the DVZ jump. Moreover, would one consider not the static potential, but 
instead the case when uj = k^ (with J± still zero), the limits of (46 1 would be different for different ways of 
bringing masses to zero: 



J2 

J o 

1.2 

A 'll 



J2 

k 2 
II 



if first Mq — > (coincides with the massless QED case) 



if first Mi -> 



(49) 



if first M = Mi 



Lorentz eigenvalues 

The other possible choice of the matrix / in Q is the Lorentzian unit matrix, which means in the case under 
consideration that one has to diagonalize the kinetic operator 



-k 2 



M 2 







— toklt 

j 2 - Ml 









Ml) 



(50) 



5 It sounds as follows. Suppose one considers the static potential at small values of mass parameters. Then, if the coefficient in 
front of fe^ in the denominator of the propagator (=an eigenvalue) is not going to zero with mass, nothing drastical happens. If, 
however, it goes to zero, one needs some further inspection of the situation. 
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Diagonalizing this matrix leads to the Lorentz eigenvalues, the result reads: 

\ L _ = \ (A L - yjAl + 4(M 2 cj 2 - Mlk\ MlMl)) (51) 
A+ = \ (A L + y/Al + 4(M 2 cj 2 Mlk\ - MlMl)) (52) 

with 

A L = iv 2 - kf\ - M 2 - Ml (53) 
and all other d — 2 eigenvalues are equal to 

\L 2 ,,. 

There are again the same two kinds of dispersion laws, 

,.2 

Ml 



Af = lo 2 - hi- Ml (54) 



uj 2 = kl + Ml (55) 



and 



since the determinant 



det K£ = det r? w det ^ = - det (57) 
Now one can again analyze the eigenvalues for physical effects: 

tachyons: Since the dispersion laws are the same, the tachyon also emerge under the same conditions as in 
the Euclidean case. 

superluminal: Similarly, the condition for superluminals to emerge are the same, 
ghosts: The derivatives -j^ on mass shell for the Lorentz eigenvalues are 



and 



9X - - M « (59) 



du 2 k 2 (M 2 -Ml) + M$ 



v ll 

Again, the ghost content of the theory may change only at Mq = (but does not change at this point) or when 
Mq — Ml < 0. The second condition is different for the Lorentz and Euclidean eigenvalues, while the first one 
is the same. Moreover, the value of is the same in both cases provided Mi — (the counterpart of mi = 
condition in the gravity theory). 

DVZ-jump: The derivatives on mass shell are now 

"hi 



L 



dk 2 



-1 (60) 



and 



d\ ± ( 



dk 2 " \ I 



(61) 



Therefore, again both at M = and Mi = there is the DVZ jump. 

Now, the lesson is that the Lorentz and Euclidean eigenvalues give the same dispersion laws and, therefore, 
the same superluminal and tachyon conditions. Moreover, at least, at the case under consideration they gives 
rise to the same DVZ-jump condition. The only difference is in the ghost content conditions. However, even 
these latter are same provided Mi — in the vector theory case, or mi = in the gravity case. 

Note that the Lorentz eigenvalues are often simpler to use, especially in the Lorentz-invariant theories (Mq = 
Mi in formulas above). Indeed, in the latter case the eigenvalues becomes functions of only the combination 
— uj 2 + fc^ and can be calculated at the rest frame (where fc|| =0). This simplifies calculations much, and if the 
ultimate results coincide, one would prefer to use exactly the Lorentz eigenvalues. 
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Appendix II. Ghosts and tachyons 



In this appendix we very briefly comment on the terminology, used in s.5 of the main text. 

Ghosts. A typical example of ghost emerges in the theory of a vector field with the Lagrangian 

- (d,A„) 2 = -Al + A - k 2 (A 2 - A 2 ) (62) 

Clearly, Aq has a "wrong" sign in front of the kinetic term and thus energy is unbounded from below. This 
means that there can be problems in constructing a full set of normalized states, or - if the theory is adequately 
modified at the strong-field regime - with making the answers independent of this kind of modification. Such 
problems are typical for ghosts and one can naturally wish to see when they can arise. At the same time, there 
is nothing bad seen in the spectrum of the theory, if we define it with the help of the Lorenz-invariant metric 
II'- {d'tu ~ ^lnv)A 1 ' — provides the same spectrum A = — lu 2 + k 2 for all components A v . If one wants to trace 
this type of ghosts already at the level of spectral study, one should better use Euclidean eigenvalues with 77^ 
substituted by 5^, as suggested in [14]. 



However, example (62 1, though standard, is not fully representative. Ghost appears here due to the vector 
nature of the field A^, but this does not mean that all ghosts should emerge for this reason only. However, 
if the rais on- d'etre is different, switching from Lorentzian to Euclidean eigenvalues does not help, as we also 
saw in [T3] in analysis of the scalar ghosts. Actually, what matters in a complicated theory is not a particular 
criterium used to identify ghosts, what matters is how normal particles turn into ghosts and/or back to normal 
as one moves around in the moduli space A4, and to see this many different criteria can be used. In s.5 we 



mentioned the simple criterium (12 1, where in the case of scalars one can use both Lorentzian and Euclidean 



eigenvalues, as we do in the present paper and in [13] respectively, so that one can compare the results. As well 



as we can judge, criterium ( 12 1 is also in accord with [5] 



Though this has no direct relation to content of the present paper, in this Appendix we also remind briefly 
what is bad (or good) about ghosts and tachyons. Usually one is afraid of ghosts for the three reasons: 

(i) they can grow in time, 

(ii) they can have negative norms and thus violate perturbative unitarity, 

(iii) they interact badly with normal particles. 

Actually, these reasons seem different and not obligatory related to each other. 

The first reason (i) is pure classical: if one begins with the Lagrangian cuf) 2 — V(4>) with a positive potential </) 
and then change the sign of a from positive to negative, then one immediately and in accordance with criterium 



(12) obtains solutions with imaginary frequencies, which either grow or decrease in time, instead of oscillating. 

It deserves noting that within the standard framework of QFT perturbation theory one deals only with 
decreasing or oscillating solutions. Indeed, if one has a free relativistic particle with the action 

(-w 2 + e 2 (k)^\c/)(k)\ 2 dujdk (63) 

one got used to define the Feynman ("casual") propagator as follows: 

z iut dw 1 f ( 1 1 \ iuJt _ Oity^ 1 + 0{-t)e~ l£ ^ t 



Lo 2 -e 2 (k)-iO 2e(k) J \ u - e{k) - iO w + e(k) + iO J 2e(k) 



which is interpreted as a particle with the dispersion rule to = +e(k) > propagating forward in time and 
antiparticle with lo = —e(k) < propagating backwards in time. Technically, the integration contour is closed 
by adding an infinitely remote semicircle in the upper half-plane (with Im lo > 0) for t > and in the lower 
half-plane (with Im lo < 0) for t < 0. Accordingly contributing are different items, with poles lying in the upper 
and lower half-planes respectively. This very fact implies that the propagator can not grow with time at t > 
and can not grow backwards in time at t < 0: exponents are never positive. One should also take into account 
different orientations of closed contour in two cases, and factor 2iri is included into the definition of integral. 
If we now consider a more general action 

(- f(k)uj 2 + g(k)^\(j)(k)\ 2 dujdk (65) 



16 



where f(k) and g(k) can become negative at some values of k (see Fig|i]), then the same propagator becomes 
more involved, but exponential growth does never occur: 



/ 



0(0' 



+fl(-t)e 



f(k)u 2 -g(k)-iO 



g(t)e-Vl g/ - f(£) l+fl(-t)eVl a/f(S) l 

2y/\fg(k)\ 

g(Oe-Vl g/ - f(g) l+fl(t)e-V g/f(E) l 
2v/|/ff(fe)| 



9(t)e-"V 9/fm +g ^_^ e ity/ g/f(k) 



when 



when 



when 



when 



/(fc) > 
g(k) > 

/(*) < o 
9$) > 

/(fc) > 
g(k) < 

/(fc) < 
«?(fc) < 



(66) 



Thus, we see that the choice of retarded and advanced Green functions is defined by the relative signs in front 
of different components of kinetic terms with respect to auxiliary term ie||</>(fc)|| 2 , added to the action. Hence, 
this choice depends on the definition of the norm for the field. In this paper we assume that these norms are 
defined by Minkowski metric i]^, namely for the gravity field \\h\\ 2 = h iiV h a p'rf- a 'rfP = . In [2] we 

used instead a Euclidean norm = Y)~}— h^ y . The disadvantage of Lorentzian norm is that by using it 

we allow existence of fields with negative norms "by hands", what does not happen with the Euclidean choice, 
which therefore can be better for analysis of the physical properties of the theory per se - as we claimed in [14 
- already for this simple reason. However, in Lorentz-invariant theories (like that of gauge-violating massive 
vectors) the Lorentzian choice looks more "natural" (preserving Lorentz invariance), while Lorentz- violation 
is introduced as a deformation of Lorentz-invariant model and thus does not allow an abrupt switch from 
Lorentzian to Euclidean norm. Since in this paper we are more concentrated on formal approaches to the study 
of eigenvalue "bundle" over the moduli space of theories than on the physical properties of massive gravity, 
we perform all analysis in terms of Lorentzian norms and eigenvalues. For analysis of Euclidean norms and 
eigenvalues - which we think is more physically justified - see |14j . 



Another typical example of the problem (ii) is well familiar from conformal field theory: Virasoro descen- 
dants of some "good" states can easily happen to possess negative norms (and one tries to get rid of them in 
construction of "unitary" CFT models). In fact here, like in the previous example of negative norms for the 
zero-component of a vector field, the problem arises only when one imposes a symmetry requirement on the 
metric in the space of fields. In CFT example one requires that the Virasoro operators are Hermitian w.r.t. the 
scalar product, induced by the norm. However, when symmetries arc explicitly violated, already at the classical 
level - that of the Lagrangian - there is no need to impose any type of symmetry requirements on the norms and 
measures (which define quantization rules), one can always use positively definite norms on entire field space, 
and in the limit where symmetries are restored and ghosts decouple the two procedures - with invariant and 
non-invariant norms - are actually equivalent, though sometime it can be somewhat tedious to demonstrate. 
Thus it is not so simple to decide what really happens, either unitarity is broken by existence of negative-norm 
states or, more probably, an anomaly occurs: the symmetry and excitation spectrum of emerging theory is 
different from the one that one could naively expect. 

The really serious problem can be the third one, (iii). The simplest example is again that of ^o-component of 
the vector field. In order to have all norms positive in this theory one needs to define vacuum state as annihilated 
by annihilation operators a|0 >= for all the spatial components A of the vector field, and by creation operator 
ao|0 >= for Aq. This means that one actually looses a possibility to define a Lorentz-invariant vacuum. This 
could still be tolerated, but the situation becomes dramatically worse when one tries to switch on interactions. 
To the best of our knowledge so far no consistent perturbation theory was developed for interacting ghosts and 
normal particles (though one can easily believe that there are problem-free non-perturbative theories of this 
kind: say when globally energy is positively defined while one begins with a perturbation theory around a local 
maximum) . 

Tachyons. This term is unfortunately used in total contradiction with its literal meaning. Lexically, 
"tachyon" referred to superluminal propagation, but today we have to use the word " superluminal" for such 
particles, because "tachyon is used to mean something else. Actually, tachyon occurs when vacuum is pertur- 
batively unstable, then it can start to decay independently at casually-disconnected points and this can look 
like a propagation of a superluminal excitation, but physically the reason is obvious and very different from 
real super luminals. From the spectral point of view, what is now called tachyon is the pole in the propagator, 
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«1 M " N 

Figure 5: A typical example of function f(k 2 ) in Lorentz-violating dispersion relation A = —f(k 2 )u} 2 + k 2 + m 2 = 0. For most 
values of space momenta we have just an ordinary relativistic particle, while in some region in k space it becomes superluminal, 
then carries an instantaneous interaction (just like Newton-Coulomb- Yukawa potential) - this happens at points M and TV — and 
between M and N it behaves like ghost, i.e. describes advanced rather than retarded interaction. 



occurring at vanishing frequency u) = 0. The archetypical example is the case f(k) = 1, g(k) = k 2 — m 2 with 
the "wrong" sign in front of the mass term. Then, at small t 

e" lkx d 3 kduj m 2 . , , . . 

J 2 -k 2 + m 2 ^o -- [J ^ mr) + N ^ mr)] (67) 

where Jfc and are the Bessel and Neumann cylindric functions correspondingly, and the prime means the 
derivative w.r.t. the argument. Therefore, in this case the propagator behaves as e" nr jr at large distances 
(small k, there is no pole), and is singular, 1/r at small distances (large k, there is a pole). The indication of a 
tachyon is non-decaying correlation at infinity. Note that simultaneously at large distances the time correlation 
exponentially decays and at small distance does not. This means that at large distances the time and spatial 
coordinates interchanged, and the causality is violated (correlations do not fall outside the light cone). 



Appendix III. Analysis of a model characteristic equation 

In order to illustrate the eigenvalue behavior it is instructive to examine a model characteristic equation which 
is quadratic but not quartic in A: 

C(A) = A 2 + (2k 2 + a) A + (/3k 2 + 7) = (68) 

In this case, one can explicitly solve all equations, and we use this to illustrate the way the information can be 
extracted from plots (which are equally available beyond the quadratic case) . The two eigenvalues are 

A ± = -( 2fc2 +^)±^ J D = (2fc 2 + a) 2 -4( / 3fc 2 + 7 ) = (2fc 2 + a-/5) 2 + (2a/3-/3 2 - 4 7 ) (69) 



Their fc 2 -derivatives are equal to 



so that 

d\+ d\- D — (2k 2 + a — f3) 2 (2a(3 - (3 2 - 4 7 ) 



dk 2 dk 2 ' D D ^ 
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Figure 6: Eigenvalues j69| l as functions of k 2 at different values of parameters a, /3 and 7. a). The case when 2a/3 — /3 2 — 47 = 0, 
discriminant D is a full square and X± become linear functions of k 2 , or X± = — ( 2fe +n)±|2fc +a — /3]^ ^ q ^ exact. The horizontal 
line is at Xca = t>). Resolution of the crossing singularity at 2af3 — /3 2 — 47 > 0, when X± are real at all values of k 2 . 

Asymptotically at k 2 — > ±00 eigenvalues tend to Aoo = — lim fc 2_ 100 ^k'^+a = ( or *° mnm ty). Punctured lines show the same 

cross as in Fig. a, the role of the resolution (deformation) parameter is played by 7. c). Resolution of the crossing singularity 
at 2a/3 — (3 2 — 47 > 0, when X± fail to be real- valued at some k 2 . d). The 3d plot of X± as function of k 2 and 7. The saddle 
structure is clearly seen. It degenerates to Fig|2a in the special case of (1756 when 2a/3 — f3 2 — 47 = 4(d — 1)B 2 and is never negative. 



while the resultant in the numerator of (20 1 is 



resultant^ 



2k 2 + a f3k 2 + 7 



2 









= -(2a(3-f3 2 - 4 7 ) 



(72) 



Denominator of ( 20 ) is simply D 



This implies that 



dx 



in two cases: either when k = ±00 and discriminant \D\ = 00 or when 2a/3 



(3 — 47 = and D is a full square, so that \± become linear functions of k . This property is nicely illustrated 



by the plots in Fig 6 In this way one can extract information from eq.(20l and from plots, what is especially 



useful in the realistic case (32 1, when Ci{\) has degree 4 and analytical approach is less straightforward. 

Explicitly in the linearized Lorentz-violating gravity ([!} the characteristic polynomial for the scalar modes 
is equal to 



C 4 (A) = A 4 + 



((d-3)(- 



ml 



m\ 



2m\ + (d- l)m 



A 3 + 



(73) 



Y- (d- 2){-uj 2 + k 2 ) 2 - ((d-3)(mjj -m\- m%) - (d - l)m 2 ^)uj 2 + ((d- 3)(mg - m 2 - m% +m%) - 2{d - 2)m|)fc 2 

2 2 o2 2,o2 2, 4 , / j 2 2 2 2 2 2 4\\\2 

—m in 1 — 2m m,2 + 2m 1 m 2 + m 2 + (d — l)(m m 3 — m^^rn^ — m 2 m 3 — ra 4 ) J A — 



(d - 2) [(mj - m 2 )u> 4 - 2{m\ - m\ + m\- m%)Lu 2 k 2 - (m 2 + m 2 - mj^fc 4 ] - 

— ((d — 3)(r?iQm 2 + m^m 2 — m 2 m 2 ) + (d — l)(r?iQm| — m\m 2 — m|))w 2 + 
+ - 3)(mgmf + m\m\ - m\m\ - m\m\ + mfm 2 + m|) - 2(d - 2){m\ + m|)m|)fc 2 
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-2rnQ'infni2 — vri^m^ + m\m\ + (d — l)(m\ + m\){m^ ) m^ — m|) — (d — l)mf jn|rn|J A+ 



+ — 2) ymQm^Ld — 2(mQjTt 2 — m m s — m-^m^ + m 4 )ui k — m 1 (m 2 — m s )k J — 

-m\((d — 3)mQm2 + (d — l)(mQm 2 — m|)^toi 2 + m\ Ud — 3)(m,Qm 2 — m^m^ + mf) — 2(c/ — 2)m 2 m. 4 ! ) fc 2 

2 2 4 i u ,\ 2 2/ 2 2 4\' 

-rajirijmj + (a — l)m 1 m 2 (mo m 3 — m 4) 



At A: = 0, i.e. in the rest frame, the characteristic polynomial (73) factorizes: 

C 4 (A)| £=0 = (A + ^ 2 -m 2 )(A-TO 2 ) 



• |a 2 + \{m\ - ma+ (d- !) m 3 - ( d ~ 2 )^ 2 ) ~ TO o( TO 2 + (d- 2)cj 2 ) + (d - l)(mo7 
For ([2]) the last bracket turns into 

C 2 (A) = A 2 + \{dB - 2A - (d - 2)w 2 ) + ((d - 2) (A - B)uj 2 + A(A - dB) 



it)} 



with the two roots given by the last row in ( 29 ) , 



A. 



A 



dB-(d- 2)oj 2 ± y/(d - 2) 2 (B - uj 2 ) 2 + 4(d - 1)B 2 



(74) 



(75) 



(76) 



Example (68 1 can be now used upon identification fc 2 = — ^^uj 2 , a = dB — 2 A, (3 = 2(B — A), 7 = A(A — dB) 
and 2aP — p 2 — 47 = 4(d — 1)B 2 > 0. In general in the rest frame one gets from the last bracket in (74 1: 



a = m\ — to 2 . + (d — 1)7713, j3 = 2m 2 ,, 7 = (d — l)(mQm§ — to 4 ) — TOqTO 2 . and 2a/3 — /3 2 — 47 = 4(d — 1)to 4 > 0. 
This implies that the singularity is resolved at ui 2 = in the single possible way, which explains the universal 
structure of Fig[IJ>. 

If instead of fc = we put u> = 0, characteristic polynomial (74(A) factorizes in a less radical way: 



a 



(A)Uo = (A-to 2 ){a 3 + ((d-3)fc 2 + m 2 -2m 2 + (d- 1)to 2 )a 2 + 



(77) 



+ ( - (d - 2)fc 4 + [(d - 3) (to 2 , - mj + to 2 ) - 2(d - 2)to 2 ] fc 2 + [to 2 - 2to 2 to 2 + (d - IXtoqTO 2 - to 2 to 2 - m 4 )] ) X+ 
+ ((d-2)(m 2 -7n 2 Jfc 4 +[(d-3)TO 2 ^TO 2 .-m 2 .)+2(d-^^ 

The roots of this equation are plotted in Figj4] Of course, in the Lorentz-invariant case ^ this ([77]) further 
reduces to (A — fc 2 — A) times (75), with —uj 2 substituted by fc 2 . 

For the full C 4 (A) in (73 1 the resultants of C 4 (A) with and 
complicated and essentially different, however, they contain two common d-independent factors 



au2 and "^2' m t ne numerator of (|20| are rather 



and 



{{m^ + TO 2 ) (TO 2 — TO 2 + TO 2 ) — mfj 

I 2 4 / 2 , 2 2\ 2i*2 2^4 / 2 2 . 22, 2 2\ 2 , 

I TO 4 W — (TOq + TO 2 — m 3 )UJ fc — TO 4 fc — (TOqTO 4 + TO 2 TO 4 + TO 3 TO 4 )CiJ + 

,/22, 22 4 r>4\T2 , / 2 2 2 , 222 6 \ 

+ (TO TO 3 + TO 2 TO 3 — TO 3 — ZTO 4 )fc + (TO TO 3 TO 4 + TO 2 TO 3 TO 4 — TO 4 ) 



(78) 



(79) 



The first factor ( 78 1 vanishes in the Lorentz invariant case ([2| , when one of the eigenvalue lines is obligatory 
horizontal. Furthermore, rcsultantA ^C(A), 9 g^2 - ) ~ fc 2 , while resultant^ fc(A), ^i^- 1 ^ ~ cj 2 , so that they 

vanish at fc = and uj — respectively - in accordance with Fig s[T]d and [2]d, which both contain one horizontal 
eigenvalue line (associated with the spatial Stueckelberg scalar). In addition, in the rest frame, at fc = 0, the 
second factor ( |79[ ) is proportional to m 4 , what corresponds to appearance of the second horizontal eigenvalue 
line in Fig[lji at m 4 = 0. As to the first factor (78 1, when it vanishes beyond the Lorentz-invariant case ([2]), it 
still signals that a horizontal line occurs. In fact, this is the line A = m\: if (78) vanishes, then C±{X) is divisible 
by (A — to 2 ) for all values of uj and fc. 
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Appendix IV. Eigenvalue bundle over the moduli space 



1 Deformation of four crosses 



The four-cross pattern, FigjTJi, describes the dependence of four scalar eigenvalues of 
k 2 — 0. Then the two horizontal lines are \ s s = —m\ and As = to 2 , while two other lines with slopes +1 
(normal particle) and — (d — 2) (ghost) are given by \ s t = —uj 2 + m 2 and \ st T — (d~2)uj 2 + m 2 — (d— 1)to|, i.e. 
they intersect the ordinate axis at m\ and m\ — (d — 1)to 2 respectively. There are two propagating (on-shell) 
modes with A = 0, one is normal, another one is ghost. 

When m\ ^ is switched on, FigJTjD, one of the four crossings is resolved: the one between X s s and X s tr- 
The intersecting eigenvalues repulse but intersection with the abscissa axis corresponds to a ghost in both cases, 
m\ > and m\ < 0: depending on the sign of to 2 the on-shell ghost comes from either lower or upper of the 
two branches. The only exception is the case of m\ = 0: then this on-shell mode simply disappears at infinity 
and the ghost is eliminated, at least at k 2 = 0. 

Switching on k 2 ^ resolves all the four crossings (even if m\ = 0). This adds one more option: that ghost 
can be eliminated not only when m 2 , = but also when m 2 = 0, exactly in the same way as abo ve. H owever, 
with increasing k 2 and m\ the patterns deviate pretty far from the four-crosses of Figfl^, see Fig 7]8 



for some 



examples. This means that even if ghosts are eliminated in the vicinity of four-cross pattern, they can re-appear 
at larger values of the mom entum k 2 . Also the case of large m| requires more careful analysis. 

As clear from Figs 7j8 the only two possibilities to have ghost-free models arise at either m\ — or at 



ml = 0. 



2 Analysis through the chain of bifurcations 

Analysis of the properties of propagating particles can be performed in a certain order, because actually there 
is a hierarchy of interesting properties. 

1) Plot A(— u> 2 ) at fixed k 2 and masses or A(fc 2 ) at fixed u> 2 and masses. Of interest are on-shell states 

A = and the slopes — J^| A _ or Jp at these points. In what follows we consider mostly the first option: 
A(-c 2 ). 

2) The signs of derivatives are actually controlled by topology of the graph A(w 2 ), especially by positions of 
the branching points: zeroes of discriminant discrim^C(A)^ where different branches merge or intersect. These 
critical points w 2 r are themselves not on-shell, but they actually define the properties of on-shell modes. They 
depend both on k 2 and masses. Of primary interest is their dependence on space momentum k 2 at fixed masses. 

3) The properties of on-shell particles change qualitatively at bifurcation points when uj 2 r merge, vanish 

or go to infinity. This is controlled by zeroes of the next-level discriminant discrim^w 2 r (fc 2 )^ . These zeroes 

fc 2 r (masses) depend only on masses and change when the masses change, i.e. when we move along the moduli 
space M.. At some points of M there can be regions in momentum space where on-shell particles are ghosts, 
and regions where they are always normal or are ghosts for all values of k. 

4) The boundaries between these regions are defined by the next-order discriminants discriml k 2 r (masses) 



One can again make an iterative study: change first some of masses, most conveniently, to 2 , and then the others 
considering higher and higher order discriminants at each step. 



3 Examples 

We give now examples of such hierarchical analysis. 

1) Some plots of the four-branch function A(— oj 2 ) are shown in Fig sj7]and[8| In Fig [7] the values of masses 
are fixed and different plots are for different values of k 2 . In Fig(8]we fix k 2 — 1 instead, but change two of the 
five masses (to 2 , and to 2 ) instead. If another mass (to 2 ) is changed we get a very different pattern, Fig|9] There 
is no problem in making many more plots of this kind. The problem is to find some reasonable way to put this 
huge collection in order. This is what above hierarchical procedure is supposed to do. 

2) From Figs{7]and[8]it is clear that the whole pattern is very well controlled by position of the branching 
points (where the tangent line becomes vertical). These branching points can be defined by pure algebraic 
means: they are zeroes of discriminant: solutions to the equation 

discrim A (c(A)) = (80) 
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Discriminant at the l.h.s. is a little too long to present here, but it is an explicit polynomial] and can be easily 
evaluated for any particular set of masses. After that its zeroes can be found numerically and they are plotted 
in the center of Figjll] as functions of k 2 for the same values of masses that were chosen in Fig|7] One can 
easily compare Figs |7] and [IT] and the moral is that essential information about the pattern in Fig (7] is actually 
contained in the far simpler and pure algebraic plot in Fig |ll| In fact, one can easily plot zeroes of the same 
discriminant as functions of masses instead of k 2 and reproduce the essential properties of Fig (^instead of Figj^J 
3) Fig _ 

Fig 



itself can be changed if we vary remaining parameters. In Fig 

shows what happens with Fig 111] ■ 



were plotted as functions of k 



12 
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positions of the branching points 



| when one of the mass parameters is 
changed. The difference between FigsfTTland Fig(12] could be systematically controlled in an algebraic way, if 
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we look at zeroes of the repeated discriminant - the crossing/merging points of the three branches in Fig 
This structure is discussed below and pictured in Fig[T0][T3] 

4) Procedure can be repeated again and again, going to higher and higher codimension in the moduli space 
M.. Thus we obtain a systematic approach to the study of bifurcations/reshufflings of eigenvalue bundle and to 
construction of phase diagrams of the theory. 



4 The bundle structure 

It is instructive to present the same in a slightly different words - and pictures,- by making more explicit the 
structure of eigenvalue bundle over the moduli space of linearized massive gravity. 

Let us fix the values of four masses, to 2 , = 4, m\ — 0, m\ — 4, m 2 — 6 and look what happens when 
we change mf from 2.64 to 3.53. The choice of masses is rather arbitrary with two exceptions: m\ is taken 
vanishing in order to look at appearance and disappearance of the on-shell ghost, and m 2 is chosen to vary in 
the vicinity of the critical value, where 

A = m\(m\ - mf,) - m 4 = (81) 

and where an interesting bifurcation occurs. In this particular case the critical value of m 2 is mo a/ to| — m 2 , = 
2^2 = 2.828427.... 

Thus we begin from mf = 2.64. Over this point of the moduli space M. there is a fiber of our "bundle", 
consisting of the 4-branched function X(u>, k). Instead of hanging such 3-dimensional fiber over the base point 
we do another thing: we hang first a 2d plot of uj 2 rit (k 2 ), which shows how the three critical points a, (3, 7 - the 
three zeroes of discriminant discrim^C(A)^ - change with the variation of momentum k 2 . After that over each 
point of this fiber - which we call discriminant fiber in what follows - we hang the 2d plot A(— uj 2 ) (the eigenvalue 



fiber), as shown in Fig 10 In this particular case of m\ = 2.64 discriminant fiber consists of a single real branch 
and only a single branching point a is seen in the eigenvalue plot. We show also the enlarged vicinity of a in 
accompanying figure, where one non-very-interesting branch is not actually seen. It is clear from this picture 
that there is a single on-shell scalar, it is located at oj 2 = and it is ghost, because the slope of the branch is 
negative at the intersection with abscissa axis. In fact, this is a very exotic excitation, being simultaneously a 
carrier of instantaneous interaction uj 2 — and a ghost since dX/duj 2 < on shell. Its characteristic dispersion 
relation is —uj 2 k 2 = it [JJ]. Actually it remains of this same kind for all values of m 2 , only the coefficient 
at the l.h.s. becomes a sophisticated function of k and can even change sign with the variation of k 2 . The 
physical implications of such " instantaneon" excitation in the spectrum remain an interesting subject for future 
investigation. 

Now we start increasing m\. At m 2 = 2.828428 . . ., i.e. at the critical value A = a new couple of branches 
shows up in discriminant fiber (they were complex at lower value^] of mf) and they get well separated soon 
enough, we show an example at m\ = 2.9. It is see that at this value of m\ one of the two new branches merges 
with the old one at fc 2 ~ 2 - and this is reflected in the properties of the eigenvalue fiber. 

At larger value of m 2 = 3.53553 . . . the other two branches merge as well, and it is quite interesting to look 
at the corresponding collection of the eigenvalue plots, shown enlarged in Fig[TT] We see that with the change 
of k 2 the on-shell scalar converts from normal particle at k 2 < 1.623759 . . . into ghost and then back into normal 



6 As a function of its coefficients discriminant of the order-4 polynomial C(A) = ^^t—o ^i^* ' s gi ven by 

-\C±ClC\ + 16C 4 C|C - 128CjC^Cf - 27C|C| - 6C 4 C7 C*|C7f + f44C* 4 C^C* 2 Cf + 144Cf CbC 2 C? + 18C 4 C 3 ClC 2 + Cf Cf Cf- 

-4Cf Cf C*o - 4C'lcf + 256Cf C$ - 192Cf C^C 3 Ci - 8OC4 C 3 CiC|C* + f8C*f CiC^Co - 27Cf Cf 

ft is a polynomial of degree 2(4— 1) = 6 in the coefficients and can be read from the celebrated Sylvester formula or represented as a 
combination of two simple diagrams, see [18j[f9|. Substitution of coefficients expression through masses, frequency and momentum 
from (|73| makes this expression rather lengthy. 

7 One can see that the two complex branches are already very close to become real by looking at eigenvalue plot at m| = 2.64 
and k 2 = 6: it is pretty clear that something is going to happen, and this behavior of the eigenvalue curves signals that discriminant 
zero is nearby - just not seen in the real section of generic complex picture. 
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particle at k 2 > 11.39337. . .. Note that the crossing of branches in discriminant fiber at k 2 « 3.85 does not 
cause any reshuffling in the eigenvalue fiber: this is because this is crossing rather than merging of branches. 

It deserves mentioning that the slightly-virtual instantaneon can actually be described analytically: at small 
values of u> 2 the corresponding eigenvalue is 

Xinstant — ~ " ~ ~t~ O(tJ^) (82) 

(d - 2)(m| - m 2 )k 4 - [(d - 3)A + 2(d - 2)m 2 m 2 4 ) k 2 + (d - l)m 2 (m 2 m 2 - mf) - m 2 m\ 



where A is given by (81). One can easily check that this simple formula provides a full description of the 
bifurcations which we are shown in the above pictures. Thus, the exactly solvable example confirms the results 
of generally applicable discriminant analysis. 



23 




24 



Figure 8: The graph of the dependence of eigenvalues on a; 2 for k 2 = 1 and different values of masses mi and mo- Masses at 
the left pictures (down from above) are equal to m 2 , = 0, m 2 = 2, m 2 . = 4, m| = 6, m| = 0; m 2 , = 10, m 2 = 2, m 2 = 4, = 6, 
m? = 0; rrij = 10, m 2 = —2, m 2 = 4, = 6, rrt| = 0. Similarly, those at the right pictures are nig = 10, m 2 = 0, m 2 = 4, 
m| = 6, m| = 0; m\ = -10, raj = -2, m| = 4, mf = 6, m| = 0; m 2 , = -10, m 2 = 2, m 2 , = 4, m\ = 6, m| = 0. Therefore, the 
two upper pictures correspond to the case when one of these masses is equal to zero. 
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Figure 9: Variation of the eigenvalue curves in Fig[7]with the change of mass 1714: m| = 3; 2.82843; 2.64575 (from the left to 
the right). The scale in the middle picture is different. 
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